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Ground state and spin-wave spectrum of cubic magnets with the Dzyaloshinskii-Moriya interac- 
tion (DMI) such as MnSi and FeGe are studied theoretically. Following interactions are taken into 
account: conventional isotropic exchange, the DMI, anisotropic exchange, magnetic dipole interac- 
tion and Zeeman energy. In the classical approximation these interactions determine the helix wave 
-vector k and critical field H c for the transition to the "ferromagnetic" spin configuration. This field 
depends on the sample form due to the demagnetization. The linear spin-wave theory is developed. 
The spin-wave spectrum depends strongly on the magnetic field. At H > H c we have quadratic 
spectrum with the gap linearly increasing with the field. Below H c the spectrum is gapless and 
strongly anisotropic. It is a result of incommensurate magnetic structure when the DMI breaks 
the total spin conservation law and umklapp processes appear connecting the spin-wave excitations 
with momenta q and q ± k with different energies. For q along k the spectrum is linear. For other 
q directions it have very complex form determined by solution of infinite set of linear equations 
connecting the states with q and q± nk where n — 1, 2, .... Restriction to n = 1 gives six equations 
which general solution remains complex. For q _L k there are two modes: one has the gap equal 
to Ak 2 y/2 where A is the spin-wave stiffness at q 3> k. The second is gapless and proportional to 
qj_. At q± ^> k all branches merge and the anisotropy of the spectrum disappears. These results 
changes insignificantly in n = 2 approximation. 

The classical energy depends on the field component along the helix axis k only. However it is 
known experimentally that rather weak perpendicular field rotates the helix and its axis is settled 
along the field. This quantum phenomenon is a consequence of the umklapps too. The spin-wave 
spectrum is unstable at infinitesimal perpendicular field. If the gap A is introduced the spectrum 
becomes stable if g/J,BH± < A In this field there are the magnetization along Hi and deformation 
of the helix. The gap appears due to cubic anisotropy and the spin-wave interaction considered in 
the Hartree-Fock approximation. 

Peculiar properties of the ESR and neutron scattering in the helical magnets are considered and 
possibilities of corresponding experimental studies are discussed. 

PACS numbers: 61.12.Bt;75.30.K 



I. INTRODUCTION 

Itinerant cubic magnets MnSi, FeGe, FeSiCo, etc have attracted a lot of attention last years due to their specific 
electronic and magnetic properties. The former are characterized by closeness to quantum phase transition, which is 
achieved at high pressure (seei^ and references therein) . Corresponding theory was developed in£. The last is related 
to the P2i3 symmetry which allows the Dzyaloshinskii-Moriya interaction (DMI) responsible for magnetic helix 
structure. In cubic crystals the DMI fixes the the sense of the helix (right or left-handed spiral) but can't determine 
its direction. It is stabilized by very weak anisotropic exchange interaction (AEI) onljii anc£. The helix structure is 
very sensitive to external magnetic field. Two field-induced transitions were observed^ and£ and references therein. 
In low field it is transition to the state with the helix axis along the field. Then, there is the second transition 
to the "ferromagnetic" state. These transitions were described by the mean- field theorjui which contains a lot of 
phcnomcnological parameters. It should be noted also that critical properties of MnSi near transition temperature 
T c are very unusual too. Corresponding experimental studies and its theoretical explanation may be found in RefsiiSii 
andi£. 

In this paper we present theoretical description of the low-temperature properties of the cubic magnets with the 
DMI. We microscopically evaluate the ground-state energy and spin-wave spectrum. In particular we demonstrate 
that the critical fields mentioned above are related to the parameters of the spin-wave spectrum in agreement with 
existing experimental data. 

The outline of this paper is follows. In Sec. II the theoretical model is formulated. Along with the exchange 
interaction , the DMI and AEI considered ir£ we include also the magnetic dipole interaction and the Zeeman 
energy. Following^ we introduce all three spin components in each lattice point of the helix and demonstrate that in 
"ferromagnetic" state the perpendicular spin components responsible for the spin- waves remain rotating. Sec. Ill is 
devoted to consideration of the ground state energy. The critical field for the transition to the " ferromagnetic " state 
is determined with the demagnetization corrections. In Sec. IV the spin- wave Hamiltonian is considered. It consists of 
two parts: direct and umklapp. The first part is a conventional spin-wave Hamiltonian. Its diagonalization determines 
the spin- wave spectrum. The second one describes transitions between spin- waves with momenta q and q ± k and 
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different energies. These transitions are a result of the incommensurate helical structure and low symmetry of the 
DMI. This umklapp interaction is considered in Sec.V. It is shown that the umklapps created by the field perpendicular 
to the helix wave-vector k are a cause of its turn to the direction of the field. At zero field the umklapps bring to 
strong anisotropy of the spin-wave spectrum: excitations with momentum along and perpendicular to k have different 
energies. Spin configuration in weak perpendicular field is discussed in Sec. VI. Origin of the spin-wave gap A is 
studied in Sec. VII. There are two contributions to A 2 : the spin- wave interaction considered in the Hartree-Fock 
approximation and the cubic anisotropy. If the last contribution is negative and sufficiently strong instead of the 
helical magnetic order the chiral spin liquid state is realized. Possibilities of the ESR and neutron scattering studies 
of the dynamical phenomena investigated in the paper is discussed in Sec. VIII. The principal results are summarized 
and discussed in Sec. IX. Some details of calculations are presented in Appendixes A-D. 



II. MODEL 



There are following principal interactions in cubic magnets without center of inversion: i) conventional exchange 
interaction, ii) the DMI, iii)the AEI, iv)magnetic dipole interaction and Zeeman energy. We investigate spin con- 
figuration in magnetic field and the spin-wave spectrum using these interactions. Then we discuss role of the cubic 
anisotropy in formation of the spin-wave gap. 

There are four magnetic ions in the unit sell of the magnets under study. We will consider here the total spin of 
the unit sell as a basic magnetic entity completely neglecting internal movements of the cell spins and corresponding 
optical spin-wave branches. Interested in spin dynamics, we can use the magnetic density approximation^ or the 
effective spin-lattice model with total spin of the unit cell S(R) and usual commutation relations 

[5„(R),5 fl (R)] = ie a/37 5 7 (R) (1) 

The DMI destroys ferromagnetism and is responsible for the long periodical spin-density wavei^j^. For each spin we 
have right-handed local orthogonal coordinate frame with basic unit vectors (r, t)r and £r which for helical structure 
are given byi 3 - 

Cr = c sin a + (a cos k • R + b sin k • R) cos a, 

?}r = —a sin kR + b cos kR, (2) 
£r = c cos a — (a cos k • R + b sin k • R) sin a. 
where a x b — c, b x c — a and c x a = b. 

Corresponding projections of the spin operators have well known form 
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here 5 is a parameter connected to the cell magnetization by M — g^sS/a 3 where [1b > and g ~ 2. For MnSi 
magnetic moment per spin 0.4^b is strongly reduced in comparison with the paramagnetic state value 1.4/is and the 
unit-cell spin 5 = 0.2 x 4 = 0.8l . 

For the above-mentioned principal interactions in the R space we have 
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It is convenient to represent these expressions in q space 

Hex = 2 ^2 J i S i ' S -i> 
Hdm = i L> q q ■ [S q x S_ q ] , 

H ae = \ ]T FWuS^S\, (5) 

v — x,y,z 

^ = y£(S q -$)(S_ q -$), 

H z = A rl/2 H-S . 

where N is total number of the cells, S q = N~ 1 ' 2 SRexp(iq • R), q — q/q, uj a = An^g^s) 2 /vq is the characteristic 
energy of the dipole interaction and vq is the unit cell volume and we omitted the isotropic part of the dipolc 
interaction which has the same symmetry as the exchange one. The DMI and AEI are of the first and second order 
of the spin-orbit interaction respectively. So we have J 3> D 3> F. Value of the dipole energy loq will be given in the 
last Sec. 



III. CLASSICAL ENERGY 



Replacing Sr by S'Cr we get the classical ground-state energy per unit cell in the following form 

S 2 - S 2 F 
Ed = - ~2~(Jb sin 2 a + J k cos 2 a) - S 2 D (k ■ [a x b]) cos 2 a H — [/c 2 (a 2 + b 2 x ) 



(6) 

+ ky (a 2 + by) + k 2 {a 2 + b 2 )} cos 2 a + Sh\\ sin a + '^—uj N cc sin 2 a, 



S 2 



here hii — g/J-H^ , Hu is the field component along the c-axis and N cc is the corresponding component of the demagne- 
tization tenso*l&. It is important to note that the classical energy depends on the field along c axis only. Meanwhile 
experiment shows that the system is very sensitive to weak perpendicular field&ji and&. We will explain below the 
nature of this quantum phenomenon. 

We are interested by small k only (ka <C 1) when Jk = Jq — and = Dq. As we will see below A is the 
spin- wave stiffness constant. At Hu = this ground state energy was studied in [4]. From Eq.(6) for components of 
the vector k we get 



from this expression we obtain 



SF 

Ak v + —K(al + bl) = SD [a x b] v . (7) 



Ak 2 + — I(k) = SD {k ■ [a x &]), (8) 



where I = ^2 k 2 (a 2 + b 2 ) is a cubic invariant. Inserting k v from Eq.(7) into Eq.(6) we get 

f {S 3 D 2 S 3 F 2 I 2 S 2 FI\ , „, . S 2 N CC . 2 
= \ 2A + 8 + I a + 11 Sm " + ~~2^ U ' a ' ^ 

where the F 2 term may be omitted. There are two possibilities^. For negative F the classical energy (6) is minimal 
for k along one of the cubic diagonals and / = 2fc 2 /3. If F > the minimum is for k along the cubic edge and 1 = 0. 
In both cases the spins rotate in the plane perpendicular to the vector k and 

_ SD [a x b] 

A + SFI/2' [W) 

For Do > and Do < we have the right and left-handed helix respectively. The classical energy depends on the 
field projection onto the vector k (c-axis) only and from Eq.(6) for we obtain 

)H\\/Hc, H\\ < H c ,,,, 
since = < " 11 (11) 

— 1, Hu > H r , 
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where the critical field is given by 

SF 

gfi B H c = h c = Ak 2 + —I + Su> N cc , (12) 

here the last term is a result of the demagnetization and the intrinsic critical field is determined by H^ nt = 
H c — 4ttN cc M, where magnetization in the high-field "ferromagnetic" state M is determined as Att/ibM = Suio 
and gfiH^ nt = h 1 ™ 1 = Ak 2 + SFI/2 ~ Ak 2 . For H > H c we have ferromagnetic spin configuration, but according to 
Eqs. (2) and (3) rotation persists in the perpendicular, spin-wave components of the spin density. 



IV. SPIN-WAVES 



We are interested in the long-wave spin-waves with q < k. At q 3> k one can neglect all interactions exept the 
exchange one and the spin- wave spectrum must be the sane as in ferromagnets: e q = Aq 2 . In the low-q region the 
hierarchy of the interactions becomes very important. The ferromagnetic exchange J is the strongest one. The DMI 
and AEI are result of the weak spin-orbit coupling A. The former and the last are of order A J and A 2 J respectively. 
According to Eq.(lO) k <~ A. 

We are dealing with incommensurate spin structure where umklapp processes are important which mix excitations 
with momenta q and q± k and different energies.lt is useful to represent the unit vectors given by Eq.(2) as 

Cr = csina+ (Ae ik ' R + A*e~ ik ' R ) cosa, 

77r = ^Ae ^k ■ R - iA*e~ iUR , (13) 
£ R = ccosa- (Ae tk ' R + A*e^ k - R )sina 

where A = (a - ib)/2, A ■ A = A* ■ A* = 0, A • A* = 1/2, A x c = -iA, A* x c = iA*, and A x A* = ic/2. Using 
this representation we obtain S q = S° c + S^A + S A A* where 

5 q = S* q sin a + 5 q cos a, 

Sq =S$_ k era a- Sj_ k sin a + iSjJ_ k , (14) 
S q+k sma-iSl + 

In the DMI interference terms appear proportional to S^S A k and S^S A +k . They are responsible for umklapps 

with q — > q ± k. At the same time the processes with q± 2k appear in the S A S A and S A S A terms . Due to algebra 
of the vectors c, A and A* these double umklapps are in Hae and Hd only. Besides as we see below the perpendicular 
field initiate processes with q ± k too. 

Here we evaluate the direct contribution to the spin-wave energy. The umklapps will be considered in the next 
section. 

The exchange energy does not contain the umklapps and has the form 

Hex =~ 1( J q sin2 a + J i< k cos2 a ) S i S -ci + J ^ S q S \ + ( J q cos2 a + J aM sin2 «)S|Sl q 

+ i(S*Siq - S*S%)JV q , k cosa + (J q - J q ,k)(S* q Si q + S q Si q ) sin a cos a - i(S%S^ - S* q 5^ q )7V q . k sina], 

(15) 

where J q , k = (J q+k + J q _ k )/2 ~ J - A(q 2 + k 2 )/S and JV q , k = (J q+k - J q _ k )/2 ~ -2Aqk/S. 
Direct part of Hum is given by 



H DM = -D (k[a x b} ) ^[S q Si q cos 2 a + S*S1 9 + S* q Sl q sin 2 a - (S q Sl q + S* q S£ q ) sin a cos a] 
+ *A>I>[fi x b])[(SiS\ - ^Si q ) cosa + (S^ - SlS\) sin a], 
where we neglected q-dependence of the DM interaction. For the anisotropic exchange we have 

H AE = ^ J2 {(c,g,) 2 (5 q sin a + S|cosa)(5l q sin a + Si^ cos a) + (g ' + + ^ 

ci : iy—x,y : z 

[(S^cosa — S| sina)(S_ q cosa — S_ sina) + S^S 7 1 } — 2iq u k v (a 2 u + cos a — S| sin a)S!*}, 



(16) 



(17) 
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In the case of the dipole interaction we obtain 
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^{(<7 • c) 2 (5 q sma + S , q cosa)(S , £ q sina + 5l q cosa) 



2|q,k| 



:[(S^ cosa — S' q sina)(S'£ q cosa — S£ q sina!) + S^Sl^}, 



(18) 



where q\ = q 2 + q 2 and |q, k|~ 2 = (|q + k|~ 2 + |q — k|~ 2 )/2. In this expression we have taken into account that the 
vector k is perpendicular to the ab plane. 

Replacing in Eqs.(15)-(18) 5 q by N^-^SSqfl we obtain a sum of terms proportional to Sq which cancels corresponding 
contribution from the Zeeman energy. In the sum Hex + Hdm + ^ae the terms proportional to sin a and cos a 
cancel too due to Eq.(8), if one neglects small terms of order of Fk ■ q. 

We consider now linear spin-wave theory. From Eqs. (3), and (15)-(18) follows 



H sw = i E i a t a i + -rr( a q a -q + a q a ± q )L 



where 



Ak 2 



£ q = Aq 2 + (— + 



SFI S 



1 ' -)cos a H — — [g c cos a 



Q i 

- — f-pr(l + sin 2 a)] — S(h c sin a + h\\) sin a, 

|q,k| 



(19) 



(20) 



where hn = gixsHn and h c is determined by Eq. (12). Using Eq.(ll) one can show that expression in the last brackets 



is zero if Hn < H c . For £> q we have 



Ak 2 SFI 2 Su 2 
(— + — ) cos a +—(lc 



|q.k| 



r ) cos a, 



The spin- wave energy is determined as e q = (E 2 — B 2 ) 1 / 2 . For Hn > H c we have 



e q = Aq 2 + h\\ - h c + 



At Hn < H c the energy has the form 



Aq 2 + 
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(21) 



(22) 



(23) 



and we have the gapless excitation with linear dispersion at q <C k. In these equations the factor q\ /|q, k| 2 is singular 
if q — > ±k and should be replaced by q\ /2 or N±/2 = (N aa + Nbb)/2 if q_L = 0. This replacement becomes evident 
if we put q = k + q^ . 



V. UMKLAPP INTERACTION 



There are following umklapp contributions to the Hamiltonian: i) Interaction of the perpendicular magnetic field 
with S 1 * spin component. ii)Intcrference of the longitudinal (c) and transverse (A and A*) spin components in the 
DM and dipole interactions. All these contributions mix excitations with q and q ± k. The dipole interaction contains 
also the (A A) and (A* A*) terms which mix waves with q and q± 2k. We will show below that this contribution 
may be neglected. The umklapps in the anisotropic exchange is neglected too. 

The umklapp contribution to the Hamiltonian consists from terms odd and even in operators a and a + . It is 
convenient to write out general expression for it and then extract the bilinear part which contribute to the spin-wave 
spectrum and linear one considered in the next section. So we have 

H^ = h • A(Sl k sin a — iS^ k + a q+k a q cos a ) 

— h • A*(5 k sin a + iSQ + a q _ k a q cos a), 

H U = £{[-2Do(q • A) + cj (c ■ q)(q ■ A)]S q S A q 

+ [2,Do (q ■ A*) + u) (c ■ q)(q ■ A*)]S q S A q }, 
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where H v = + H$P , and S c , S A and S* A * are given by Eq. (14). In this expression we included also terms 

linear in S_ k . They are not contribute to the spin- wave energy and will be considered in next section. 

The term S5q t o in the expression for leads to demagnetization of the perpendicular field if H\\ 7^ which may 
appear due to low symmetry of the sample. In this case — h • A is replaced by 

P = -(h- A + Slj N cA sin a), (25) 

where N c a — (N ca — iN c b)/2 is corresponding off-diagonal component of the tensor N a p^. It is an unusual feature 
which has to be clarified. The tensor N a p is diagonal if the coordinates are along the principal axes or the ellipsoidal 
sample. If k is not along one of these axes the off diagonal components of N a p appear. The 5 q ,±k terms in S A ' A do 
not contribute due to the condition q • A = 0. As a result for the bilinear spin- wave contribution we obtain 

H$w = ^{[f cosa + i?_ q _k(l - sin a) 
— i? q (l + sin a)]a q " +k a q + [P* cos a 

+ iT q (l - sina) - i?*_ k (l + sina)]a+_ k a q (26) 
+ i2q[a q a_q_ic(l — sin a) — at q Oq +k (l + sin a)] 
+ -R- q [a± q aq_ k (l - sina) - a q a_ q+k (l + sina)]} 

where 

= ^Y^[-2D ( q ■ A) +u>o(c- q){q ■ A)). (27) 

Here again i? q contains off diagonal component of N a p. From these equations we see that the umklapp interaction 
contributes to the spin-wave energy at Hu < H c only and do not affect excitations with q along k. At the same time 
from Eqs. (15)-(18) for the direct interaction follows that at Hu > H c the zero-point vibrations disappear. 
For consideration of the umklapp interaction we will use equations of motion 

wGab(^) + G [h , A ],b{u) =< [A B \ >> (28) 

where H = H^, + Hgfy and the Green functions are determined by 

G a ,b(") = ~* / < [M*)> B (°)} >=-<A,B> u (29) 



where the expression in the right hand side coincides with the conventional determination of the generalized suscep- 
tibility. 



We will use below the following notations 



"q'°q M a q±k' a q (30) 

F a+ ,a+ = F 11 F a+ „ ,a+ = F ±l- 



-qTk' 



Neglecting the umklapp interaction we have 



G q (w) = —5—^; F q ( W )=-^-i T , (31) 



w - e q or - e q 



where Eq, Bq and e q are given by Eqs. (20)- (23). If Hu > H c we have F = 

Using Eq.(28) we obtain infinite set of equations which contain along with conventional Green functions G and F the 
functions G±„k and F± n k where n = 1,2.... Below we will restrict mainly by n = 1 approximation. Its validity will 
be discussed later. As a result we obtain six linear equations for functions (30) (see Appendix A). In general form 
their solution is very complex. So we will consider two limiting cases: q — > and q > k. 

The q — > case. Corresponding equations are analyzed in Appendix B. Neglecting the dipole interaction we obtain 
the final result for weak perpendicular field in the form 

(Ak 2 /2) cos 2 a + uj 



Gq( W ) 



2 -el + (h 2 /2)cos 4 a' 

(32) 

_ . . (Afc 2 /2)cos 2 a V ' 

Fq(CJ) 



e 2 +(h 2 ± /2) cos 4 a' 
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where h"j_ — (gfiB) 2 (H 2 + H 2 ). We neglected here the small terms proportional to h\ in numerators as well as the 
contribution of the slightly splitted modes with the energies e ~ ek (See Appendix B). The functions G± and F± 
are of order of hi and are evaluated in Appendix B. In these expressions the square of the spin- wave energy in the 
transverse field at q = is given by Eq.(B5) 

e q (H x ) = A 2 - (hl/2) cos 4 a. (33) 

For the gapless spin- waves Cq(Hi) is negative for any Hi and magnetic subsystem is unstable. It means that the 
helix axis k must turn and stands along the field. For finite gap A, the spin-wave spectrum remains stable up to 
hi ~ A and then there is the first order transition to the parallel state. In the intermediate case, when Hn ^ there 
is rather complex behavior governed by the equation of state derived in the next Sec. Similar situation takes place 
in conventional antiferromagnets: If the field is along the sublattice magnetization the first order spin- flop transition 
occurs at H = A where A is the spin-wave gap. For the inclined field there is the rotation of the sublattices along 
with the first order transition^. We consider origin of the gap below. 

The rotation of the helix axis was observed in FeGe— and MnSA and&. In both cases Hi is much less than H c . 
For MnSi the transition to parallel state is at Hi c ~ 0.1T and H c ~ 0.6T. It is the reason of the approximation 
with n = 1. Any additional n produces small correction of order of H\jH 2 -C 1. 

The q^k case. As above we consider the first-order umklapps connecting q and q± k. We consider q 1 k case 
only where the umklapp interaction the dipole interaction is absent (see Eq.(27). We put also H = 0. Then we 
discussed briefly the role of small qn. We demonstrate also that n — 2 approximation changes does not change 
qualitatively the n = 1 results. 

If q _L k we have two interacting modes with q and q ± k which energies according to Eq. (23) are given by 

e q = Aa(<£ + fc 2 ) 1/2 
e*±k = A[tf x + k 2 )tf x + 2k 2 )] 1 '*. 

The excitation e q is the gapless and linear in q at q <C k. The gap of the double degenerated mode is equal to Ak 2 \[2. 
The umklapp interaction lifts this degeneracy (see Appendix C) and we have three different modes. One of them 
remains non-renormalized and coincides with the gapped mode in Eq.(34) and two other are given by 

e ± = A[k 2 + Ml + ql/k 2 ± (fc 2 + 8q 2 ± + \lq\jk 2 + 8q 6 x / 'fc 4 )] 1 / 2 , (35) 

where the plus mode has the same gap as in Eq. (34) but other dependence on qi . 

The non-renormalized mode is irrelevant as corresponding expressions for the Green functions given by Eqs.(C5) 
do not have poles at this energy. It is evident that for other directions of q the Green functions must have all three 
poles. The gapless mode is strongly renormalized and at small qi we have quadratic dispersion with e_ = Aq x /2. 
Asymptotically at q k all three modes tends to Aq 2 and the umklapps become unimportant. The umklapp 
renormalization is seen clearly in Fig. 1 where the full and dashed lines correspond to minus and plus modes in 
Eq(35) and the dash-dotted line to gapless mode with q || k respectively. 

If q || k the umklapp interaction disappears (see Eq. (27)). General expression for arbitrary directed q is very 
complex. But it becomes very simple for small qi u <C k. In this case from Eq. (C3) we get 

e cl = A[qi/2 + k 2 q 2 ] 1 ^. (36) 

Cancellation of the q x term is the most striking feature of this expression. So we have to discus validity of our n = 1 
approximation. Corresponding equations of motion are given by Eq. (Al) where at H = and q = q^ and the 
amplitude of the umklapps i? q = — S'Do(q ■ A). So at small qi <C k and lu -C Ak 2 we have an expansion in powers of 
q\ as G± and F± are proportional to 1/ek- So the n = 1 approximation gives correct result for q 2 ^_ term at small u 
and above cancellation of q\ survives in higher approximations. It shown in Appendix C for n — 2. However the q\ 
remains with slightly changed coefficient and instead Eq.(36) we obtain 

e q = A[3qi/8 + k 2 q 2 ]^ 2 (37) 

It is follow from the above argumentation that further approximations can not change coefficients in this expression. 
For qi 3> k we have 1/qi expansion and the umklapp renormalization disappears. Unfortunately at qi ~ k and 
lu ~ Ak 2 the results of n — 1 approximation are qualitatively correct only. It is demonstrated in the end of Appendix 
C. 
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VI. PERPENDICULAR SUSCEPTIBILITY AND HELIX DEFORMATION. 

We have demonstrated that the terms linear in operators a and a + cancel in the diagonal part of the interaction 
given by Eqs. (15)-(18) due to equilibrium conditions for the classical energy. However there are not similar conditions 
for the terms with a±^ and a^ k in the interaction (24) with the perpendicular field and we have 

+ + (38, 

This expression may be considered as perturbation which gives rise addition contribution in the energy as well as some 
non-zero parts of the average values of the transverse spin components < Sff 1 > proportional to h ^ . But we will use 
more general method, which leads in the first approximation to the same results. We consider operators a±k and a± k 
as c-numbers. Corresponding c-number terms appear in Eqs. (19) and (26) and we obtain additional contribution to 
the ground-state energy 

Ei = P^Jj2[a+(l + s) - o_ k (l - a)] + [Pc - R Q {1 + s)]a+a Q 

[Pc + R (l - s)]a£a- k + i? [aofl-k(l - s) - aja+(l + s)] + h.c. (39) 
+ £k(|a k | 2 + | a _ k | 2 ) + B k (a k a_ k + a+a\) + E Q \a \ 2 + B (a 2 Q + a +2 )/2, 

where s = sin a, c = cos a, i?o = (S/2)loqN c a cos a and in the classical limit a + — a*. 

The energy E\ is minimal at the following obvious equilibrium conditions dEi/dd^± k = which lead again in 
n = 1 approximation to six linear equations. Corresponding matrix coincides with M given by Eq.(B2) at u> — and 
we get equation 

M(0)Z = V, (40) 

where column Z — (a , aj, a_ k , at, a k , a^ k ) and 

V = (0, 0, -P* y/S/2{l - s), ~P*^/Sj2{l + a), P^J r SJ2{l + s),P-y/S/2(l - a)). (41) 

In the case when dipole contribution Rq = solution of these equations is presented in Appendix B. Using above- 
mentioned equilibrium conditions one can show 

Ei = | ^1 [K + a£)(l + a) - (a_ k + ol k )(l - *)] + c.c, (42) 
and for the field depended part of the ground-state energy we have 

F = J_H sj^tf 

Ak 2 2Afc 2 (l + cos 2 a)[A 2 - (/i 2 L /2)cos 4 a]' k ' 

where the first term is the Zeeman energy for H < He- 

The helix axis has to rotate at g/j,BH± ~ A. According to experimental data&^ the rotation was observed at 
H±_ <C H c . So we can put cos 2 a — 1 and obtain 



Sh 



,,, 2 

E ~ Ak 2 



9 sin 2 ^ 
cos" * 



i(l-h 2 x sin 2 ^/(2A 2 ) 



(44) 



where VP is the angle between the helix axis k and the field. For small h± this energy is minimal at \& = i.e the 
field along k. The real direction of the helix axis is determined by competition between the magnetic energy and 
anisotropic exchange given by Eq. (9). Eq.(44) is not valid at h ~ A as it was derived in the linear approximation 
when the amplitudes a and a + given by (B7) should be small. 

The perpendicular field deforms the helix structure producing the higher harmonics. In our n = 1 approximation 
there is the second harmonic only. The static contribution to the transverse spin components has now the form 

S R = iAe lkR (^e jkfl + S\e~ ik - R ] + c.c, (45) 
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where iS± k = y / S'/2(a±k — a ^ k ) an< ^ we P u ^ a = 0. From these expressions using Eqs. (B7) we obtain the transverse 
uniform magnetization^ 



$U = - TTJT^xf ( 46 ) 



5hj.A 2 
2Afc 2 (A 2 -h\/2) 

The second harmonic of the helix structure induced by the perpendicular field is given by 

9A 2 

S * = -2M^Wm [MA ' ^ + ^ ■ ^ 1 (4?) 
As we have mentioned above corresponding second-order Bragg reflections were observed 6 and 8 . 

VII. CUBIC ANISOTROPY AND THE GAP PROBLEM 

We postulated above existence of the spin- wave gap. Now we demonstrate that there are at least two contributions 
to the gap: cubic anisotropy and interaction between spin-waves considered in the Hartree-Fock (HF) approximation. 
The former may have arbitrary sign and the last is positive. So different contributions to the gap may compete. 
Changing the sign and strength of the cubic anisotropy by pressure one can get the quantum phase transition from 
the ordered to the spin-liquid state. This phenomenon was observed in MnSi^f^ and^i. 

Complete expression for the gap at arbitrary H\\ is very complex. So we present here results for H = and H > H c . 
We begin with the cubic anisotropy. As is well known the single ion cubic anisotropy has the form 

V = Kj2iSt, R + S£ R + Si R ). (48) 

R 

Using Eqs. (3), (13) and (14) for zero magnetic field (a — 0) we obtain 

V = 6S*K \A„\ 4 + 6S 3 K {2|A t/ | 2 c 2 a+a k + [|A„| 2 c 2 -2|^| 4 ]( aka _ k + a + a ; k )}. (49) 

v=x,y,z k,u=x,y,z 

Neglecting in Eq. (9) the F 2 term for the classical energy we have 

E " = ^T E kl{al + bl) + ^ y, + (5°) 

u—x.y,z v=x,y,z 

The spin- wave terms in Eq.(49) give additional contribution to i^q and in Eq.(19) which are given by 

SE, 1 = 12S 3 KJ2Z 2 M V \ 2 ; «? q = 12S 3 i^(c 2 K| 2 -2K| 4 ), (51) 
and for the anisotropic contribution to the spin-wave gap we obtain 

A 2 c = 3 -S*Kh c £ (al + bir, (52) 

v=x,y 1 z 

where h c — Eq + Bo, h c , Eq and Bq are given by Eqs. (12), (20) and (21). This equation for the gap holds at q = 
when g 2 — > N cc . For small q we must replace h c by hl nt + Sw^q 2 . So at q = and q ^ the gap depends on the 
sample form and angle between q and the helix axis respectively. We see that A 2 is positive for K > only. In 
ferromagnets this sign of K corresponds to the easy directions along the cubic diagonals. 

Let us consider now the classical energy and the gap for orientations of the helix vector k along [1,1,1], [1,0,0] 
and [1,1,0] directions labeled as 1, 2, and 3 respectively. We can choose ci = (1, 1, l)/3 1 / 2 , ai = (1, -1, 0)/^ 1 / 2 
and k = (1, 1, -2)/6V 2 ; £ 2 = (1,0,0), a 2 = (0,1,0) and b 2 = (0,0,1); c 3 = (1, l,0)/2 1 / 2 , a 3 = (1, -1, 0)/2V 2 an d 
63 = (0, 0, 1). For the classical energy we have 

S 2 Fk 2 S A K 3S 4 K S 2 Fk 2 9S 4 K . , 

EcX\ = — 7. 1 7T~ ! E cl2 = — - — ; E c i 3 = — 1 — — ; {06) 

6 2 4 8 16 

From these expressions we see that structures (1,1,1) is realized if S 2 Fk 2 < 3S 4 K/2. Otherwise we have (1,0,0) 
structure. These expressions are derived at T — 0. However interplay between two anisotropy constants F and K 
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must determine the helix direction at all T and explain the transition from (1, 0, 0) to (1, 1, 1) structure in FeGe with 
decreasing T— . 

Corresponding expressions for the gap are given by 

A 2 cl = 4S 3 Kh c ; A 2 c2 = 6S 3 Kh c ; A 2 3 = (9/2)KS 3 h c . (54) 

We consider now contribution to A 2 arising from the spin-wave interaction considered in the HF approximation. 
It is well-known that in the Heisenberg ferromagnets there is not a gap due to the total spin conservation law. The 
DM interaction brakes this law. From Eqs. (15) and (16) we have 

Hex + H d DM = -\ ^{[J q , k - J + 2£ q (q • [a x S])](S£s! q + S«S1J + (J q - J )S q Sl q } (55) 

where we have taken into account that S q S_ q = NS(S + 1) and the odd terms with S^S^ surviving at arbitrary q 
are zero in the HF approximation. We replace also D by D q as the HF contribution is saturated by large q. Using 
now Eq.(3) for the spin- wave interaction energy we obtain 

Vj = -i ^2{V 1 [a^ +2 a 2 a^_ 1 a 3 + ~(ai - at^a^^a^] - ^(Ji - J )(a x + at-i)a^ +2Jt3 a-2a a } , (56) 

where 1, 2, 3 = qi, 2 ,3 and V\ = J^k - Jq + 2D 1 (k ■ [a x b]). 

The HF contribution to A 2 is calculated in Appendix D by the method used in22 for evaluation HF corrections to the 
spin-wave energy in 2D CuOi planes where the gap appearing due to pseudodipolar interaction between neighboring 
Cu 2+ ions. From Eqs. (52) and (D7) we obtain the final result 

a 2 = 5>j + id 2 + ^ E % ( 57 ) 

where the last sum has to be less than unity. 

In strong field H > H c we have cos a = and zero-point fluctuation disappear (B q = in Eq. (21)). As a result 
one can show that there is not the HF contribution to the gap and the classical energy and the gap are given by 

2 o (58) 

A = h\\—h c + AS 3 KY,hci + -cl(al + S 2 )]. 

For the cases listed above cubic contributions to the classical energy and the gap are given by 

F - KSA - F -K^- F KS4 - 

^cll — — 5 — i &cl1 — AO , &cl3—Z — , 

SKS 3 (59) 
Ai = — — ; A 2 = -AKS 3 ; A 3 =KS 3 . 



VIII. EPR AND NEUTRON SCATTERING 



Both phenomena the EPR and the inelastic magnetic scattering are described by the spin susceptibility. We outline 
here some of its features. For simplicity we use the n — 1 approximation. 

From Eq.(13) and definition (29) we obtain following general expression for the susceptibility 

(60) 

+ < , S_„ >uj A a Ap-\- < , iS^q >w A a Ap, 

where — < A, B > UJ is the generalized susceptibility chiAE^)- In the first three terms both spin components have 
opposite momenta and describe the direct processes . All other are the umklapp terms. We will use below the linear 
spin- wave theory and consider two limiting cases H = and H > H c . The main attention will be paid to direct part 
of the susceptibility. The umklapp part will be commented briefly in the end of this Sec. 
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Zero- field case. 

There are now two contributions to the susceptibility along the k and in ab plane. For the first using Eqs. (14) and 
(31) we have 

Xcc(q,w) = -(5/2)[G q (w)+G q (-w) + F q (w)+i!'+(a;)], (61) 

where due to non-Hermiticy of the interaction (56) we must to distinguish between F and F + . In the EPR case when 
q = we obtain 

5A 2 

xM = ~ M^ 2 -A 2 r (62) 

This susceptibility describes respond to the external ac field. Connection between external and intrinsic susceptibilities 
is given by well known equation (see for example^*) 

X a( i = xi n p - ^oxi^N^x^, (63) 

and for x™* we obtain the same equation as for \ cc with replacement h c and A 2 by corresponding intrinsic quantities. 
It should mention that there has to be very small contribution of the t\ mode (more precisely e+ mode) . It appears if 
in Eqs. (B9-10) we retain terms proportional to (Eq — Bq)\R\ 2 = A 2 \R\ 2 /(Eo + Bo). Moreover in n = 2 approximation 
the €2 mode has to appear too. But we neglect both these contributions. 

For 9 ^ k the susceptibility is strongly anisotropic in q space. For q || k the umklapp interaction is zero and from 
we obtain 

SAqj 

Xcc(g||,tj) = - a a , (64) 

911 

where e 2 ( = Aq 2 [A(q 2 + k 2 ) + Su q 2 }. 

With increasing of q^ tree different modes appear. However at qu = there are two modes (see Appendix C) given 
by Eqs. (35) and using Eqs. (C5) we get 



SAq 2 

X cc (q±,w) = 



e 2 + -e 2 + A(2ql + k 2 ) e 2 - e 2 _ - A(2q 2 ± + k 2 ) 



(65) 



where e\ is given by the second line in Eq. (34). In the n — 2 approximation additional modes appear and so on. 
However amplitudes of these modes decrease with n. 

Consider now the second two terms in Eq. (60). By the same way as above we get 

where we put a Q a^ + b a b/3 = 5 a p — c a cp and 00,6,3 — cLpba = e a /3jC 7 . Perpendicular and chiral susceptibilities are 
determined as 

X-Lq = Xq+kH +Xq-k(^); C q (w) = Xq+k (w) - Xq-k (w) , (67) 

where 

X Q M = | [F Q («) + F+ («) - G Q H - G Q (- W )] . (68) 

Functions Imx± and ImC determine parts of the neutron scattering cross section independent on the neutron polar- 
ization Po and proportion to it respectively (See for example^). The chiral contribution to the cross section appears 
due to the DM interaction and is q odd in agreement with general theory 2 ^. 

In the EPR case when q — * the singular term q^_/|q, k| 2 in Eqs. (20) and (21) is equal to iVi/2 = 7V aa + jV&& 
and q 2 = N cc respectively. At the same time the umklapp interaction \R\ = Su>o(N 2 c + A r2 c ) 1 / 2 /4 is small and we can 
neglect it in the dispersion Eq.(B4). As a result using Eqs. (20), (21) and (31) we obtain 

X a p M = ~ C a Cf3) z » ( 69 ) 
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where 



ei = [2{Ak 2 f + Scu Ak 2 + (Slu ) 2 N cc N ± /2} 1 / 2 . (70) 

and C = 0. If we take into account the umklapp interaction the resonances appear at lu = A and £2k and so on. 
If q ^ as above we consider two cases: q || k and q 1 k. In the first case we have 



X-L 



Z 4 



c 



(71) 



where Z± = A(k ± q\\ ) 2 + Slu and e 2 ± = A{k ± q\\) 2 [A{k ±q\\) 2 + Ak 2 + Soj ] . 
For ql k using (C5) we obtain 



S 



4(e 



(72) 



where e± are given by Eq. (35), Z± = ±[2A 2 q 2 ± (2k 2 + qj)} ± (e| - e?) and C* = 0. 

We discuss now briefly off-diagonal terms in Eq.(60). They are a result of the umklapp interaction and has to be 
less than that considered above. At q = the cA and AA terms are proportional to off-diagonal components N ca ^ 
of the demagnetization tensor in the first and the second degree respectively. Possibility of their experimental study is 
a special problem which is beyond scope of this paper. At q ^ these terms are very complex and we do not analyse 
them here too. 

"Ferromagnetic" state (H > H c ) In the longitudinal cc channel there are the two-spin excitations only which are 
beyond our consideration. By the same way as above for the transverse susceptibility we obtain 



X±af3 = -~(5 a p 



C a Cp) 



1 



£ q+k 



£q-k 



1 



£ q+k 



+ 



1 



£q-k 



(73) 



where the second term represents the chiral contribution to the susceptibility and the spin-wave energy is given by 
Eq.(22). At q = the last term in (22) contains the demagnetization and = Ak 2 +Sio {N aa +N bb )/2+g^ B {H~H c x)- 
In this case the imaginary part of the chirality is cj-even and the chiral contribution to the static susceptibility is equal 
to zero in agreement with the general theory^. Experimentally this contribution could be measured using circularly 
polarized AC field. For q^ the neutron scattering is maximal at q = ±k in contrast to conventional ferromagnets. 



IX. SUMMARY AND DISCUSSION 



We begin with a short survey of the main results presented above. We have used the following interactions: 
conventional isotropic exchange, the DMI, anisotropic exchange, magnetic dipole interaction and Zeeman energy. In 
the classical approximation these interactions determine the helix form and the critical field H c of the transition to 
the "ferromagnetic" spin configuration. This field depends on the sample form due the demagnetization. 

The linear spin-wave theory was developed. It was shown that the spin-wave spectrum depends strongly on the 
magnetic field. At H > H c we have quadratic spectrum with the gap linearly increasing with the field [see Eq.(22)]. 
Below H c the spectrum is gapless. It is strongly anisotropic due to incommensurate helical structure and low symmetry 
of the DMI. As a result umklapp processes appear which connect the spin-waves with momenta q and q ± k and 
different energies. For arbitrary q the energy has very complex form. It is determined as a solution of infinite set 
of linear equations connecting the states with q and q ± nk where n — 1,2, .... Restriction to n = 1 leads to six 
equations which general solution remains complex. For qlk there are two modes given by Eq.(35). One has the gap 
equal to Ak 2 where A is the spin- wave stiffness at q ^> k. The second mode is the gapless with quadratic dispersion 
at small q. Both gapless branches (parallel and perpendicular to k) and the gapped one are shown in Fig.l. At q 3> k 
all branches merge and the anisotropy of the spectrum disappears. The n = 2 approximation do not changes these 
results qualitatively. For q <C k at H = the spectrum has simple form given by Eq.(36). 

The classical energy depends on the field component along the helix axis k only. However it was shown experimen- 
tally that rather weak perpendicular field H± <C H c rotates the helix and its axis is settled along the field^. 8 . This 
quantum phenomenon is a consequence of the spiral spin structure when the angle between spin and the field depends 
on the lattice point. As a result we get the umklapps again and the spin- wave spectrum in the gapless case is unstable 
at infinitesimal perpendicular field. So we must introduce the gap A "by hands" and the spectrum becomes stable if 
g/j,BH± < A [see Eq.(33)]. The magnetization along Hj^ appears and the helix becomes deformed (see Eqs.(46) and 
(47) respectively) . The second harmonic with the wave- vector 2k was observed in& and£. 
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We considered two ways of the gap origin: cubic anisotropy and interaction between spin waves. The first contri- 
bution is proportional to the strength K of the cubic anisotropy. The second appears is a result of breaking of the 
total spin conservation law by the DMI. It is positive and disappears in the "ferromagnetic state" at H > H c due 
to absence of the zero-point vibrations. But in this region there is the field-induced gap determined by Eq.(22). In 
ferromagnets the sign of K determines the direction of the easy axis. In our case the helical structure is stable if 
A 2 given by Eq.(57) is positive. Otherwise we have the chiral spin liquid. So there is a question : If change sign of 
A 2 is a reason of the transition to the disordered state in MnSi at high pressure? This problem demands further 
experimental and theoretical study. 

The helix structure leads to peculiar features of the ESR and neutron scattering. In conventional magnets the ESR 
frequency is equal to the spin-wave gap. In the helical systems due to the umklapps along with this frequency there 
are more higher resonant excitations corresponding to the spin- waves with q = nk where n = 1, 2, .... AtH^O the 
chiral channel appears. It may be observed using circular radiation. In neutron scattering below H c three modes can 
be studied. The chiral channel exists at H = and q ^ due to axial-vector nature of the DMI. 

Detailed experimental work was done in the case of MnSi compound only. We now compare some of known 
experimental results obtained at ambient pressure with the predictions of our theory and discuss possibilities of the 
further experimental studies. The principal parameters are: lattice spacing a = 4.558A, T c ~ 29K, k ~ 0.035A -1 , 
saturated magnetization M = OA^s/a 3 ^ 0,016T, 4nM = 0.20T, critical field H c = 0.5 -r- 0.6T 7 and spin-wave 
stiffness A ~ 52meUA 2 — . From these data and Eq.(12) we obtain H c ~ Ak 2 /(g/iB) — 0.55T. This value coincides 
with experimentally observed critical field. For more precise comparison one must measure all parameters including 
the demagnetization N cc using single sample. 

To the best of my knowledge the ESR in MnSi was studied only io2&. Several resonances were observed but 
only one was studied qualitatively as a function of the magnetic field. Its frequency in zero field is equal to 0.93T. 
Using Eq.(70) and taking into account that Sloq = ingfiBM we obtain 0.85T. The agreement is within the error 
bars. According to Eq.(22) at H = H c the frequency is close to 0.6T ~ H c and then increases linearly with H. 
This behavior coincides with results of Ref3i. Further experimental studies are necessary. First of all it is essential 
to observe the gap A and its dependence on the perpendicular field, which accordingly with Eq.(33) should be 
A = [A 2 — (l/2)(g/isiJj_) 2 cos 4 a] 1 / 2 . The observation of the umklapp resonances £2k and 63k would be important 
too. 

The theory developed in this paper explains partly some experimental findings. It was developed in the linear spin- 
wave approximation. The spin-wave interaction was used for evaluation contribution to the gap only. We also did not 
evaluate some experimentally observed quantities such as specific heat and the site magnetization. Corresponding 
results will be published elsewhere. 
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APPENDIX A: EQUATIONS OF MOTION 

Using Eqs.(25), (27) and (28) in the n = 1 approximation for functions (29) we obtain 

(w - E)G -BF- [Pc + fl_ q (l - s) - ,R q _ k (l + + (fl-q + #q-k)(l + s)F_ 

-[P*c + R*_ q _ k (l - a) - J£(l + s)]G + - (i? q + i?I q _ k )(l - s)F+ = 1, 
BG + (w + E)F + (i?_ q + i? q _ k )(l - s)G_ + [Pc + ,R q _k(l - s ) - i? -q( 1 + S )] F - 

-(Rl + iT q _ k )(l + s)G_ + [P*c + i?*(l -s)- iT q _ k (l + s)]F+ = 0, 
-\P*c + i£ q (l - s) - i?*_ k (l + a)]G - {R\ + i?*_ k )(l - s)F + (w - E_)G_ - S_F_ = 0, ' ' ' 

-(iC q + i?*- k )(l + b)G + [P*c + i? q _ k (l -a)~ fl- q (l + s)]F + + (w + EJ)F- = 0, 

-[Pc + i?_ q _ k (l - s) - i? q (l + s)]G + (i? q + i?- q - k )(l + s)F + (w - E+)G+ - B+F+ = 0, 
(i? q + J2_ q _ k )(l - s)G + [Pc + Rq(l - s) - fi_ q _ k (l + s)]F + B+G+ + (w + E+)F+ = 0, 

where s = sin a, c = cos a, E{B) = E(B)q, E(B)± — E(B)q±k and functions E and B are determined by Eqs. (18) 
and (19). 
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In matrix form these equations are given by 

M(q,w)=J (A2) 

where / = (1, 0, 0, 0, 0, 0). Determinant of the matrix M is even function of lo and has the following general form 

Det[M] = (to 2 e 2 )( W 2 e 2 + )(c 2 _ £ 2_ } (A3) 

We obtain three renormalized energies as roots of this equation. As Det[M] is the denominator of the Green 
functions they have three poles. However at q _L k two initial spin- wave energies are equal and only two renormalized 
branches are physically relevant. Hence the Green functions must have two poles. That it is the case is demonstrated 
below in Appendixes B and C. There are five renormalized branches in the n — 2 approximation (see Appendix C) 
and so on. 

General expressions for e q0 and e 2 ± are very complex. We consider below two main cases: i. Small q <C k and ii. 
H = 0, q~k. 



APPENDIX B: THE q -> CASE 



From Eqs. (20) and (21) we have 



= Aq 2 + (Ak 2 cos 2 a)/2 + (Su /2)q 2 c cos 2 a, 



B q = (l/2){Ak 2 + Sujofc) cos 2 a, 



E x = Ak 2 (2 + cos 2 a)/2 + (Sto /2)[q 2 cos 2 a + (q 2 ± /2)(l + sin 2 a)], B x = (l/2)[Ak 2 + (Sto /2)(q 2 c - q\/2)] cos 2 a, 

(Bl) 

where Ex = -E q ±k and B\ = B q ±k- From Eq. (27) we get i? q = (Slo /2)(c ■ q)(q ■ A) cosa and i? q ±k = 0. If q = we 
must replace q 2 , q\ and (c • q)(q ■ A) by N cc , N aa + and (N ca — iN c b)/2 respectively. The matrix M has now the 
form 



lo - E n 



E n 



-Pc-R(l-s) R(l + s) -P*c + R*(l + s) -R*(l-s) \ 

R(l-s) Pc-R(l + s) -R*{l + s) P*c + R*(l-s) 

P*c-R*(l-s) -R*(l-s) lo-Ex -B x 

-R*(l + s) P*c-R*(l + s) Bx to + E 1 

-Pc + R(l + s) R(l + s) lu-Ex 

\ R(l-s) Pc + R(l-s) Si 

For the determinant of this matrix using Program Mathematica 5 we obtain 

Det[M] = {lo 2 - e\){{uj 2 - e\){uj 2 - e 2 ) - 4|P| 2 (£;o£i + - \P\ 2 cos 2 a + lo 2 ) cos 2 a 






LO + El 



(B2) 



8|i?| 2 (£ - B )(£i+5i) 

+ A{PR* + P*R)[lu 2 + (E Q - B a )(E 1 - Si) - 2|P| 2 cos 2 a] sin a cos a - 8|i?| 2 (£ - B a )(E 1 - B^sin 2 a}+ 
4[(PR* + P*R) 2 lu 2 - 4|P| 2 |i?| 2 e 2 ] sin 2 acos 2 a, 

(B3) 

where 6q = E 2 — Bq — > A 2 . Here the last term may be represented as {|P 2 ||i?| 2 (w 2 — e 2 ) + [(N c ±_ ■ Pi) 2 — 
N^jJP^lo 2 } sin a cos 2 a. The lo 2 term breaks rotational invariance in the ab plane and e 2 is not a root of this 
equation. But this term is very small and we neglect it. 

We assume that P, R and A are small in comparison with Ak 2 which is main parameter of the theory. In this 
approximation we have 



Det[M] = (lo 2 - ej) 



(lo 2 - 4)(lu 2 - A 2 ) - 4|P| 2 (£ £i + BoBt + lo 2 ) cos 2 a - S\R\ 



, Ei (1 + sin 2 a) + B\ cos 2 a 
Eq + Bq 



(B4) 

. As complete expressions for resonant modes and the Green functions in this approximation remain very complex 
we consider two limiting cases: R = and H = 0. 

The R = case . As we see below the Green functions do not contain pole at lo 2 = e\. For two other energies we 
have 



2 h\ (Ak 2 + Slo N cc ) cos 4 a 
* ~ 2(Ak 2 +SWVi/2 cos 2 a) 



A 2 



h\ cos 4 a 



4 + h 



(Ak 2 + Slo N cc ) cos 2 a 



2Ak 2 + Slo N±_ cos 2 a 



2 2 h 2 , (2 + cos 2 a) 2 
cos a ~ e x H — cos a, 



(B5) 
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where e\ is given by Eq.(67) and in right hand side we neglected the dipole contribution. 
Corresponding expressions for the Green functions are given by 

G(w) = Z~\{u + E )(lo 2 - e 2 ) + (hl/2)(E 1 - u)], 

F(lu) = -Z^BoIlu 2 -el- (h 2 L /2)B 1 cos 2 a], 
G-(w) = Z- 1 (h+/2)[uj 2 -el + E E 1 + B B 1 - (hi/2) cos 2 a] cos a, 
F-(w) = -Z- 1 (h+/2)[B Q E 1 +B 1 E + uj(B 1 -Bo)] cos a, 

where Z = {uu 2 — e 2 + ){u 2 — e 2 ,), h± = (h a ± ihb)/2 and one gets G+ and F + replacing h + by h- in the expressions for 
G_ and F- . Neglecting the dipole interaction we obtain (32) 

Solution of Eq.(39) for the spin deviations frozen in the perpendicular field is given by 

_ + _ fS (A • h)(l + cos 2 a + sin a) A 2 
Ok -(a,,) - -y 2 Afc 2 (l + cos 2 a)[A 2 - (h 2 ± /2) cos 4 a}' 

A* • h) (1 + cos 2 a - sin a) A 2 
2 Ak 2 (l + cos 2 a)[A 2 - (h 2 ± /2) cos 4 a] 

S h\{EQ — Bo) sin a cos a 



_ + _ o ^ ■iij^i^cos (t-smap / R? x 
a -k-ia-kJ - V IT TT27; ; —2 Nr A 2 7T2 /o\ — _1> ^ 



" ! W 2 2Afc 2 [A 2 -(/7,i/2)cos 4 a] 
The H = case. From Eq.(B4) wc have now 

Det[M] = (lu 2 - e\){u 2 - e 2 )(co 2 - e 2 + ), (B8) 

where 



c 



1 - 



{Ak 2 + Saj N cc ){2Ak 2 + Slu N ± ) 
A 2 (Sa Jo ) 2 (N 2 a + N 2 b ) 



(B9) 



2 _ 2 , " W^u; yiv ca T ^" eb 

e + — e l + "" 



(Afc 2 + Suj N cc )(2Ak 2 + Slu N ± ) ' 
Corresponding expressions for the Green functions are given by 

G(lo) = Z-\uj + Eo){u? - el) + A\R\ 2 {E 1 + Bi)], 
F{u) = -Z-^Boiu 2 - e 2 ) + 4\R\ 2 (E 1 + B,)}, 
G-(w) = R*(lu + E - B )(lu + Ei+ Bi), 
F_(w) = -R*{uj + E - S ) (Si + B 1 - u), 



(BIO) 



where R = (5 , w /4)(7V ca - iN cb ), G+ = -(R/R*)G- and F + = -(R/R*)F_. Note that now it is addition sign minus 
in expressions for G + and F + in comparison with the R — case. 

APPENDIX C: THE q > k CASE 

. It is convenient now to use dimcnsionless variables X = uu/{Ak 2 ), Y = |<fj_|/A;, V = q\\/k and U = —(q a — i%)/k. 
In these variables we have M = (Ak 2 ) 6 m and for m we get 

(X-W 2 -\ -\ U U* \ 

2 

u* 



[/ 

V o u o o i x + w 2 + 2v + §,/ 

(CI) 

where W 2 = Y 2 + V 2 . 

General expression for Det[m] is very complex. For V -C 1 we have 

Dct[m] = (X 2 - e 2 )(X 2 - e 2 + )(X 2 - e 2 _) 

e\ = 1 + AY 2 + Y A ± (1 + 8Y 2 + 17F 4 + 8F 6 ) 1/2 + 4V 2 



_ 1 


u 





u* 





X + VF 2 + i 





u 





£/* 





X - W 2 + 2V - | 


_ 1 








[/* 


i 2 
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X + VF 2 - 2F + | 

















X - W 2 - 2V - | 


_ l 


u 








3 ^ " 


rT4^ 2 + 
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where e 2 = (Y 2 + 1)(Y 2 + 2) is the energy at ± k in the dimensionless units. As above at Z = the factor 
X 2 — e 2 cancels in the expressions for Green Functions. 



For small Y and V we have 



X 2 _= Y — + V 2 ; 



2 + 4Y 2 



Asymptotic expressions for e± at 7 > 1 are given by 

e± = Y 2 ± V2Y. 

The Green functions in dimensionless units are given by 

G(X) = Z~ X \{X + Y 2 + 1/2)(A 2 - e 2 ) + 2Y 2 (Y 2 + 3/2 - X)]; 
F(X) = -Z-\X 2 -e 2 + 2Y 2 )/2; 
G_(X) = -Z-\Y a + iY b )[-(X + Y 2 + 1/2)(A + Y 2 + 3/2) + 2Y 2 + 1/4]; 
F-0 = Z-\Y a + iY b ){2X - l)/2, 



(C3) 



(C4) 



(C5) 



where e 2 = (Y 2 + l)(Y 2 + 2), Z = (X 2 - e 2 _)(X 2 - e\) and G+(F+) = [G-(F-)](Y a - iY b )/{Y a + %Y h ). 

To illustrate that n = 1 results are at least qualitatively correct we consider now the n = 2 approximation. In this 
case the matrix M = (Ak 2 ) 10 m and for m we have 
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Y 9 
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(C6) 



where A T = A =p y 2 . The determinant of m is given by 

Dei[m] = 1600A 2 - 1760A 4 + 564A 6 - 44A 8 + A 10 + (7760A 2 - 4868A 4 + 524A 6 - 15A 8 )y 2 
+ (-600 + 9279A 2 - 2462A 4 + 167A 6 - 3A 8 )y 4 + (-585 + 3686A 2 - 597A 4 + 22A 6 )y 6 
+ (-519 + 866A 2 - 107A 4 + 3A 6 )y 8 + (-528 + 128A 2 - 7A 4 )y 10 + (-144 + 24A 2 - A 4 )y 



(C7) 
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We see that according to argumentation in Sec. 5 we have cancellation of the Y 2 term. There is the y 2 A 2 term only. 

The equation Det [m] = determines now five spin-wave energies generated by e 2 = y 2 (y 2 + l), e 2 . = (y 2 + l)(y 2 +2) 
and €2k = (y 2 + 4)(y 2 + 5). We consider below the three lower branches only. 

Small Y case . From Eq.(C6) for the gapless branch we have A 2 = 3y 4 /8 and the factor 1/2 in (C3) has to be 
replaced by 3/8. Higher approximations with n > 2 can not change this result. Consider now the branches connected 
to e k . For A 2 ~ 2 the first three terms in Eq.(C7) give 648% 2 - 5904^ + 9398y 4 = 0, where X = A 2 - 2. Solutions 
of this equation are given by 



A^ 



37y 2 

18 



2.iy 2 ; 



A 



i27y 2 

18 



2 + 7.iy 2 . 



(C8) 



We again obtain small changes of the numerical coefficients. 

Large Y. For large Y we can write A = WY 2 and asymptotically W 
powers of Y we obtain the dispersion equation 

8(1F-1) + 15(VF-1) 2 



1. Retaining only three main terms in 



(w - iy 



40 



Y 2 



Y 4 W 2 



0. 



Its solutions are given by 



A± = y 2 ± YV2; Ai 



+ 2- 



(C9) 



(C10) 
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APPENDIX D: THE HARTREE-FOCK GAP 

As the interaction (55) is non-Hermitian the spin-wave Hamiltonian have the formS^ 

Hsw = ^[£?qa+a q + i(B+a q a_ q + B q at q a+)], 
where _B* 7^ Corresponding corrections to Eq and B q in Eq. (19) are determined as 

S 2 Vi \ t „ I 5 2 V! \ r , / 5 2 Vi " 



6E a 



SdqScLq 



SB, 



Sa^SaZq 



5B+ 



SaZqSaq 



At q = we obtain 

5E ° = ~\ + - J q + J oK + (\v + y q - J q + Jo)/,]; 

5S = — ^ 5S o + = ^0 - 5 5}(V& + Vq - 2 J q + 2J )n q - (J q - J )/ q 
where expression for Vq is given below Eq.(55) and 

En P 4. _i_ 

n q = — ^- ; / q =< a q a_ q >=< al q aj >= - — 

Using Eqs. (20) and (21) and neglecting dipolar terms we get 

Ak 2 ^ 

Af = - — £)[(V + Vq)n q + (Vb - J q + J )/q]. 



(Dl) 



(D2) 



(D3) 



(D4) 



(D5) 



Large q ~ 1/a give principal contribution to the expression (D 3). However Eqs. (20) and (21) determine Eq and £> q 
at q <C 1/a only. Using Eq. (55) we obtain these functions for all q 



Eq = S 



Jo- 



Ak 2 



B a = 



S 



(D6) 



where using Eq.(8) we replaced k • [a x b] by Ak 2 /(SD ). 

The function i? q oc k 2 . As e q = (E^ — B^) 1 / 2 we have that n q is proportional to (Ak 2 ) 3 and may be neglected. In 
the / q term at large q we have Vb + Jo — J q — £q/S and 



4P . 



2^ q Afc 2 

^q,k <^q 



(D7) 



where the last term contribute to the gap only as ^ J q k = 0. 
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FIG. 1: Fig.l Spin-wave dispersion for different directions of the wave-vector q. (a) q is along the helix axis k (dashed line), 
(b) Gapless and gapped branches for qlk. Solid and dot-dashed lines respectively. 
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